Entanglement detection from interference fringes in atom-photon systems 
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A measurement scheme of atomic qubits pinned at given positions is studied by analyzing the 
interference pattern obtained when they emit photons spontaneously. In the case of two qubits, a 
well-known relation is revisited in which the interference visibility is equal to the concurrence of the 
state in the infinite spatial separation limit of the qubits. By taking into account the superradiant 
and subradiant effects, it is shown that a state tomography is possible when the qubit spatial 
separation is comparable to the wavelength of the atomic transition. In the case of three qubits, 
the relations between various entanglement measures and the interference visibility are studied, 
where the visibility is defined from the two-qubit case. A qualitative correspondence among these 
entanglement relations is discussed. In particular, it is shown that the interference visibility is 
directly related to the maximal bipartite negativity. 
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I. INTRODUCTION 

Much attention has been raised recently to the study 
of atom-photon systems, such as ultracold atoms and 
trapped ions, in view of quantum information and com- 
putation processing protocols [l| • Achieving these proto- 
cols with hybrid systems requires the capability of ma- 
nipulating atoms and photons with high precision as well 
as long coherence time for qubits. As an example of such 
demands, an ability to prepare highly entangled states is 
needed in the beginning of protocols. It is then neces- 
sary to verify the quality of these entangled states and 
this can be accomplished by performing a quantum state 
tomography Q. However, it is rather difficult to perform 
the full tomography, in general, and it becomes impracti- 
cal as the size of Hilbert space increases. An alternative 
to the quantum state tomography is to perform a set 
of measurements to detect the amount of entanglement 
in the state. It is the main objective of this paper to 
investigate a passive measurement scheme for entangle- 
ment detection by analyzing the interference pattern of 
photons emitted by atoms. 

Since entanglement and correlation are closely related 
to each other, a relationship between the interference vis- 
ibility and entanglement in bipartite systems has been 
anticipated in the context of wave- particle duality [H, 0] ■ 
Jacob and Bergou observed that two-qubit concurrence 
is equivalent to two-particle interference visibility in the 
case of pure states p). To be self-contained, let us con- 
sider a simplified argument for the known correspondence 
between the visibility and the concurrence by the detec- 
tion of released photons in the far-field [3, Q • Consider 
two two- level atoms in which energy levels are |e) and 
\g), respectively. Assume an initial atomic state living in 
the subspace spanned by |0) = \eg) and |1) = \ge), 



whose concurrence is C{p) = 2|poi|- The spontaneous 
emission process of a single photon (i.e., p — > \ggk)) oc- 
curs with equal probability but different phases. Ignoring 
the details of the coupling between the atoms and the ra- 
diation field, the standard Michelson's visibility is defined 
by 



V(p) 



(2) 



where 7 max = max <$< 2 7r Tr(/($)p) and 7 min = 
mino<$<2ir Tr(J($)/o) are the maximum and the mini- 
mum of interference fringes, respectively, and the pro- 
jection measurement of the phase $ is given by = 
|$)($| with |$) = (|0)+ e i$ |l))/V2. It is straightforward 
to evaluate the visibility in this model as 



V(p) = 2\p 01 \ = C(p) 



(3) 



P=(|0>,|1>) 
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and this completes the claim (fj . We emphasize that vis- 
ibility here is not the two-particle visibility defined from 
a combination of correlation functions, which is shown 
to be equivalent to the concurrence of the pure state 
spanned by \gg), \eg), \ge), |ee) 

Several generalizations addressing the case of mixed 
states, of two qudits, and of multipartite qubits have also 
been reported recently [7Hll|. Despite these progresses, 
there is still no net conclusive result on general multipar- 
tite systems. This is simply due to the fact that no gen- 
erally accepted entanglement measure and no unique def- 
inition of interference visibility exist in multipartite sys- 
tems. The former, for example, can be seen from the fact 
that none of these existing measures can be used to order 
the entangled states uniquely [l2l - fl4j . The latter prob- 
lem was addressed in Ref. [ll| where a systematic con- 
struction for the visibility and the so-called predictability 
were introduced starting from minimal requirements for 
a proper definition of these concepts. 

In this paper, we analyze the interference pattern gen- 
erated by photons spontaneously emitted by a set of two- 
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level atoms and we investigate the relation between en- 
tanglement and interference visibility. To explore various 
kinds of two-qubit and three-qubit entangled states, we 
consider identical two-level atoms pinned at given po- 
sitions and initially prepared in a superposition of the 
first excited states of the atoms. In the course of time, 
a photon is later released from the excited atom in any 
direction. In our setting, quantum statistical effects play 
no role and complications related to the atomic motion, 
such as Doppler and recoil effects, are avoided. In an 
experimental situation, this would be realized, for exam- 
ple, with atoms trapped in different potential wells in the 
Lamb-Dicke regime. 

Another motivation of this work is to clarify the mean- 
ing of the standard interference visibility ((2J defined in 
two-path interferometers when it is applied to three-path 
interferometers. In Ref. for example, it was shown 
that there is an infinite family of state functions which 
could be considered as equally good measures of the inter- 
ference strength. In the present paper, however, rather 
than studying these many other alternatives, the usual 
standard definition @ is employed to examine what one 
can learn from it. 

The content of this paper is as follows. Section II 
provides a brief summary for the known results on two 
two-level atoms interacting with the radiation field. In 
Sec. Ill the interference pattern is analyzed when the 
emitted photon is detected in the far-field and is com- 
pared with the entanglement present in the initial (pos- 
sibly mixed) state. It is confirmed that the interference 
visibility converges to the concurrence when the spatial 
separation of the qubits goes to infinity. In Sec. Ill C, a 
state tomography is shown to be possible when the sep- 
aration between the qubits is on the order of the wave- 
length of the atomic transition. In Sec. IV the previous 
analysis is extended when the initial atomic qubit state 
is pre- pared in a W-like pure state. The relationship 
is studied between the interference visibility and known 
entanglement measures. In particular, the interference 
visibility is shown to be able to detect bipartite entan- 
glement. A brief discussion is given in Sec. IV D on a 
state tomography for the W-like state. A summary and 
possible extensions of the work are stated in Sec. V. The 
Appendix contains the definitions of the entanglement 
measures used in the paper. 



II. SUPERRADIANCE AND SUBRADIANCE 
FROM RADIATIVE CORRECTIONS 

We briefly review some known results for two atoms 
interacting with the quantized radiation field. A more 
detailed account can be found in Ref. 



15] 



Consider two neutral two-level atoms which are pinned 
at given positions Xj (J = 1,2). The internal structure 
of each isolated atom consists of a unique excited state 
|e) at energy hiu e with radiative lifetime r = 1/T, which 
is separated by the energy Hujq = hcko = 2irhc/\a from 



a unique ground state \g) at energy Hujg. Throughout 
the paper, a simple scalar model for the atom-field inter- 
action is considered and any polarization effects are ne- 
glected. Fixing the origin of energy at the ground-state 
levels (i.e., at 2huj g ) the free atoms- field Hamiltonian Hq 
for this case reads 

Ha= ^ \ e i) ( e J I + H fUva k a k ' ( 4 ) 

j=l,2 k 

where a k and a\, stand for the annihilation and creation 
operators of a photon with momentum k and angular 
frequency ui = c\k\. The total Hamiltonian is H = Hq + 
V, where the coupling of atoms to the radiation field is 
given by 



v = - E d i E (&) > 

j=l,2 



(5) 



in the dipole approximation. Here the scalar radiation 
field E(x) reads 

E(x) =i^^(e ifc x a k - e- ik ^a*) , (6) 



where £ w = ^Jhuj/(2eQL 3 ) is the field strength at energy 
huj and L is the size of the quantization box used to define 
the photon modes. The dipole operator for the jth atom 



d (\ei){9i\ + \9i)(fii\) 



(7) 



Within the scalar radiation model, the dipole strength d 
relates to the radiative width of the excited state through 
r = d 2 fcp / (27rfteo) . In the following, when there is no 
possible ambiguity for labeling the internal atomic states, 
the atomic index j is omitted. 

Consider the situation where at time t = the radi- 
ation field is in its vacuum state |0), while one of the 
atomic qubits is in the excited state. In other words the 
system starts in the subspace Hp spanned by the first 
excited states without photons (i.e., \egO) and \geO)). 
Because these first excited states are coupled to the ra- 
diation field, they eventually decay to the atomic ground 
state \gg) by releasing a spontaneous photon \k). This 
process can be described by the method of the resolvent 
operator G(z) = [z — and its projection on the 

various subspaces of interest [lH ] . For example, the time- 
evolution operator restricted to the subspace Hp is ob- 
tained from a contour integral of the following projected 
resolvent operator 



Gp(z) 
V 



V- 



1 



-v, 



z-Ho-V 
|e 5 0)(e 5 0| + | ff eO)( 3 eO| 



(8) 
(9) 



where V is the projector onto the subspace Hp. To sec- 
ond order in the coupling constant, its diagonalized form 
reads 



G P {z) 



1 



1 



- fifL 



- m. 



(10) 
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where the corresponding eigenkets |±) are given by 

|±) = A={\egO) ± \geO)) , (11) 
and the complex eigenvalues Q± = w± — iT±/2 are (l6j 



(12) 

r± = T [1 ± /(for)] , (13) 
g{x) = cos{x)/x, f{x)=sm{x)/x, (14) 



where r = |xi — a?2 1 is the distance between the atoms. 

The net effect of the interaction with the radiation field 
is to lift the degeneracy between the initial atomic states 
in the subspace Hp. The physical process behind it is the 
resonant exchange of photons between the atoms which 
"glues" the atoms together and gives rise to the states 
|±) in a fashion similar to the bonding and antibonding 
states of a molecule. The new levels also acquire differ- 
ent finite lifetimes. The subspace Hp is thus irreversibly 
emptied by spontaneous emission in the course of time 
while the ground state is gradually populated. Both the 
degeneracy splitting and the lifetimes depend on the rel- 
ative distance r between the atoms. For sufficiently close 
atoms, the celebrated superradiant and subradiant be- 
haviors are recovered in which the eigenstates are given 
by |+) and |— ), respectively. The atomic dipoles are 
perfectly correlated (i.e., they oscillate in phase for the 
superradiant state and in phase opposition for the sub- 
radiant state). For |+), the radiation waves emitted by 
the two atoms interfere constructively and the system ra- 
diates more efficiently, shortening its lifetime. For |— ), 
on the other hand, the waves interfere destructively and 
the system cannot radiate, increasing its lifetime. Indeed 
one gets T + — > 2T and T_ — > in the limit k r — > fl5j . 

In the other extreme limit k^r — > oo, both lifetimes 
achieve the value obtained for a single isolated atom 
r± — > r. In this limit the atoms are no longer coupled 
by the radiation field, and the energy degeneracy in the 
subspace Hp shows up again uj± — s- u>q. 

The time evolution of the initial state decaying into 
the ground state is obtained from a contour integral of 
another projected resolvent 



Gqp(z) = Q- 



1 



H -V 



l-V. 



The matrix elements of Gqp(z) read 
ihg 



(ggk\G QP (z)\±) 



e — ifc -xi _|_ g— ifc £C2 



(15) 



(16) 



where g u = d£ u /h is the Rabi frequency at field 
angular frequency uj. One can easily check that 
(ggk\GQp(z)\— ) — > when k r — >• as expected for the 
subradiant state since it does not couple anymore to the 
radiation field in this limit. 



III. ENTANGLEMENT DETECTION: THE 
CASE OF TWO ATOMS 

A. Visibility in the infinite separation limit 

Let us now proceed to verify the result §S§ based on the 
microscopic model presented in the previous section. As 
before, consider the initial state (JIJ of atoms within the 
subspace spanned by \geO) and \egO), and parametrize 
the density matrix as 



P2 = (\egO), \geO)) : 



1 ( l + s z s x -is y \ f(egO\\ 



where s x = s sin 9 cos 4>, s y 



(17) 



ssin#sin^>, and s z 



s cos 9 are the components of the vector s with the ranges 
< 9 < it, < <j> < 2ir, and < s < 1. The concurrence 
of the state p7|) is C(p2) = ss'm9. 

In the model considered, the interference visibility ^ 
is defined by the maximum and the minimum values 
achieved by the spectral distribution of the spontaneous 
photon over all possible emission directions. This spec- 
tral distribution is proportional to the transition prob- 
ability P(p2 — > \ggk)) from the initial state (fTTj) to the 
ground state \ggk). This is obtained from the time evo- 
lution of the state p2(t) — U(t)p2U <! (t) in the long-time 
limit Tt — > oo where U(t) is the time-evolution unitary 
operator calculated from ([I7j)) as follows: 



P(P2 



\ggk))= lim Tr (\ggk)(ggk\p 2 (t)) 
rt-s-oo 



(18) 



B 4 



( w -c + )2+(r + /2)2 ( w - w _)2+(r_/2)2 



where the coefficients B± are 



l±f 

(1 ± s x )(l ± cos k-r) + — — -(s y — gs z ) sin k r 



1- 



with the shorthands / = /(fco r ) an< i g = 5(to r )j an d 
r = X\ — X2 is the relative vector connecting the two 
atoms. In Eq. (fT5|) . the spectral distribution is a weighted 
sum of two Lorentzians centered at angular frequencies 
lo± with respective widths T±. The angular frequency 
separation Aw = \uj+ — w_ | between the two peaks is 
generally small and the atoms need to be located rather 
close to each other to distinguish them. This can be seen 
from the ratio 



Aw 



g(k r) 



1 ± f(kor) 



(19) 



which is a small number when atoms are located far apart 
compared with the optical atomic transition wavelength. 

To extract the visibility of the interference pattern, the 
distribution (fl8|) is to be maximized and minimized over 
the spherical angles of the emitted photon. Rewrite (TT5)) 



P( P 2^\ggk)) = 



Tuj 



(27T) 



21,3 
ft 



+ \ C+V 2 cos(fe-r- 



(20) 
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with the notations 
l + s x 



± 



1 



( w - w+ )2+(r + /2)2 ( w - w _)2+(r_/2)2 
1+/ , i-f 



(c-c + )2+(r + /2)2 ( w - w _)2 +(r _/2)2 



l + g* 



(21) 



(22) 



and 9q = tan -1 (??/£;_). It is clear that the maxima and 
minima occur when cos(fc-r — 9q) = ±1. The far-field 
fringe visibility is thus 



e + V 2 



(23) 



It is straightforward to calculate the visibility in the in- 
finite separation limit as 



lim V(tjJ,r;p 2 ) = \ h* 2 

r—toc V 



S 2 = 

y 



= C{p 2 ) . (24) 



Therefore, the result ([3]) is valid only in the infinite sep- 
aration limit. From a practical point of view, the infinite 
separation limit is reached as soon as the two atoms are 
separated by a large distance compared to the wavelength 
Ao of the optical transition of each individual atom. 

A full characterization of the initial state p 2 is com- 
pleted as soon as the full vector s associated to p 2 is 
known, see (fTT)). Note that the angle 9q describes a shift 
of the interference pattern perpendicular to the plane lo- 
cated at equal distances from the atoms and that it con- 
verges to the angle <p in the limit r — > oo. Thus, know- 
ing the positions Xj of the atoms, the pattern shift can 
be measured in principle and 4> can be extracted from 
the data. From the visibility and the interference pat- 
tern shift, s x and s y can be obtained. However one still 
needs to extract the missing component s z from the data, 
which is unfortunately not possible in the infinite sepa- 
ration limit. Complete state tomography is studied in 
Sec. Ill C. As a last remark, the angular separation be- 
tween two consecutive fringes is inversely proportional to 
distance r. This means that the further apart atoms are 
located, the more periodic and regular the fringe pattern 
appears. Therefore, by knowing the distance between 
the two atoms, one does not need to measure the whole 
interference pattern for all angular positions. It is suf- 
ficient to record a few fringes to detect the amount of 
entanglement. 



B. Interference visibility for finite distances 

In this section, we analyze the deviation V(r, w, p 2 ) — 
C(p 2 ) of the interference visibility obtained for atoms at 
finite distances, Eq. (|23"|) , from its asymptotic value ob- 
tained for atoms far apart, that is, from the concurrence 
C(p 2 ) = s sin 9 of the two-qubit initial state. Our interest 
is in the state and distance dependency of this deviation 



and particularly in the maximal possible deviation from 
the concurrence for a given distance r between the atoms. 
The latter number can be used as a quantitative measure 
for the relationship between interference visibility and en- 
tanglement when estimating the amount of entanglement 
in the initial state from the observed visibility. 

First, from the result (|23|) it can be seen that the de- 
viation can be both positive and negative, and converges 
to zero in the limit r — > oo as was shown in (f2~4"|) . Next, 
we numerically calculate the maximal value reached by 
|V — C| as a function of r (in units of Ao) for some spe- 
cific values of the state purity s = 1.0,0.5,0.1. To this 
end, the deviation |V — C\ is maximized over the state 
parameters 9, <f> at fixed r and s. The frequency of the 
emitted photon is set at ujq. This can be achieved, for 
example, by frequency filtering in the detection process. 
Without such filtering, one would have to integrate over 
all range of frequencies in Eq. (|18|) to get the observed 
fringes. The result, depicted in Fig. [TJ shows an oscilla- 
tory decay of maxg^ |V — C| when r is increased. The 
local minima occur when k^r is an integer multiple of n, 
that is, when /(fcof) = 0. For these particular atomic 
separations, the superradiant and subradiant states have 
exactly the same decay rate r± = T but still different 
eigenfrequencies. The local maxima occur when fcor±7r/2 
is a multiple integer of 2w, that is, when g(kor) = 0. In 
this case, the superradiant and subradiant states have 
identical eigenfrequencies oj± = ujq, but achieve different 
decay rates. As one can see, the values of the local max- 
ima are insensitive to the purity of the state whereas the 
values of the local minima increase monotonically when 
the purity of the state is increased. 
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O.X 



s=l. ()(•), 0.7(o), 0.4(A), 0.1(A) 



0.6 



0.4 



0.2 



,A 




i 



4 A 



FIG. 1: Plot of maxj^ |V— C| as a function of r (in units of Ao) 
for different values of the purity s = 1.0, 0.5, 0.1. The photon 
angular frequency has been set at cjo. The local minima occur 
when f(kor) = (i.e., when r is a multiple integer of Ao/2) 
while the local maxima occur when g(kor) — (i.e., when 
r is an odd multiple integer of Ao/4). The values achieved 
by the local maxima are purity independent while the values 
achieved by the local minima increase with the purity. The 
solid line represents the analytical result for s = 0, Eq. (|25[) . 
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For s = 0, C{p2) = and we get the analytical result 



\V-C\ = 



2k r\ sm(kor)\ 
l + (k r) 2 



(25) 



It may seem counterintuitive that one gets nonzero fringe 
visibility when the purity of the initial state is s = 
(i.e., when the concurrence of the state is zero and no 
coherence seems present in the system). In this case, 
the initial state p2 describes a complete statistical mix- 
ture of the qubit states \eg) and \ge) and one can hardly 
imagine measuring interference fringes using a Young 
slit device operating with one slit shut. The reason 
is that, due to the resonant exchange of photons be- 
tween the atoms, the radiating eigenstates are not \egO) 
and \geO) but the superradiant and subradiant states 
|±) = (\egO) ± \geO))/V2. The initial state at t = is 
also equivalent to a complete statistical mixture of |±), 
and both the superradiant and subradiant states display 
well-defined relative phases between the states \eg) and 
\ge). Both states |±) lead to an interference pattern in- 
dependently in the far-field, the dark fringes of one pat- 
tern corresponding to the bright fringes of the other and 
vice verse. In addition, these states decay with different 
rates. Therefore the sum of their respective interference 
patterns does not average out except in the limit r — > oo. 
This is easily seen from Eqs. (l2~Tj) and (|2"2"j) where rj = 
and 9q = for s = and hence the visibility of the 
interference pattern for s = relates to the difference 
of the Lorentzians. For atoms that are far apart, the 
Lorentzians are essentially identical and cancel out, and 
the interference is lost. In this limit, the intuitive result 
is recovered in which the radiation takes place starting 
from the complete statistical mixture of the qubit states 
\eg) and \ge). For atoms that are close enough, however 
interference fringes do occur as the result of the overlap 
and partial blurring of two independent interference pat- 
terns, associated to |+) and |— ), with different strengths 
and fringe locations. 

As can be seen from Fig. Q] and Eq. (|23|) , large de- 
viations between V(r, w;p2) and C(p2) appear for small 
separation of the atoms, typically in the region r < Xq. 
However, the deviation can still be significant for atoms 
far apart as it only decays like 2/kor for large spatial 
separations. It is straightforward to understand why the 
argument presented in the Introduction to derive the re- 
lation between the visibility and concurrence applies to 
atoms that are far apart. Concurrence relates to entan- 
glement present in the atomic state pi spanned by the 
states | eg), \ge) irrespective of the coupling to the radi- 
ation field while visibility refers to the interference pat- 
tern obtained by coupling the atoms to the radiation field 
(i.e., it relates to the interference pattern obtained from 
the radiating states |±)). However, for atoms that are 
far apart, the resonant exchange of photons between the 
atoms is negligible and the interaction with the vacuum 
field does not lift the degeneracy of the subspace Hp. As 
a consequence, the eigenstructures to analyze radiation 
and entanglement are then essentially the same. This 



is why, for atoms that are far apart, the entanglement 
present in the initial state becomes equivalent to the vis- 
ibility of the fringes. 



C. State tomography 

In this section, it is shown that, for finite atomic sepa- 
rations, the transition probability P(p2 —> \ggk)) and the 
visibility V(r, w; pq) encode enough information to recon- 
struct the full initial state (TTT)) . This is obvious from 
Eqs. (fT8"|) and (j2li|) which can be directly used to extract 
the values s, 9, and <f> as soon as r is known. Complete 
knowledge about the initial state then provides informa- 
tion on the amount of entanglement in it. 

There arc many ways to demonstrate state tomog- 
raphy. For example, consider the simple case where 
kor > 1, but is not too large. The frequency of the 
emitted photon is set at loq as before. Expanding the 
transition probability and the phase-shift angle up to the 
first order in / and g reads 



P(p 2 -> \ggk)) 



8c 



(27rfco 



+ (s x — 2/) cosfc-r + (sy — gs z ) sinfc-r], (26) 



tan 



gs z 



VSy 



(27) 



First of all, knowing k$r calculates the values / and g. 
Second, from the probabilities at k-r = or tt in Eq. (|26[) , 
the state parameter s x is obtained. Similarly, knowing s x 
and using the values at kr = ±7r/2, the combination of 
state parameters s y — gs z can be found. Last, measuring 
the approximated phase shift 9 in Eq. (|2T|) together with 
the obtained values of s x and s y — gs z provides s y and 
thus three state parameters s x , s y , and s z . Therefore the 
full vector s is known, completing the state tomography. 

The main result of this discussion is that spontaneously 
emitted photons from the excited states (|T7|) do provide 
full information about the initial state as well as the 
amount of entanglement in it. The drawback of this to- 
mographic scheme is that it is not very efficient. Indeed, 
in the absence of any prior information on the initial 
state that could help simplify the state search, one needs 
to measure photons emitted in all directions to obtain 
the full fringe pattern and later characterize the state. 
This requires the repeated measurement of photons from 
an identically prepared state. Although the method is 
not very practical, it can be combined with other pro- 
posals. For example, other possible entanglement de- 
tection schemes were studied by many authors. Among 
them, we mention Ref. [l7T | , where an entanglement wit- 
ness was constructed without using the full interference 
fringe pattern. 
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IV. INTERFERENCE VISIBILITY IN THE 
THREE-ATOM SYSTEM 

In this section, we extend the previous analysis to the 
case of three atoms and investigate the relationship be- 
tween the visibility and the various entanglement mea- 
sures. This can be accomplished along the same lines of 
the microscopic calculation done for the two-atom sys- 
tem. As shown in the following, however, a complete 
analysis appears difficult. Two major obstacles are shown 
in the following. The full diagonalization of the resolvent 
operator does not take a simple analytical form except for 
a few limited cases, for example, when the three atoms 
are pinned at the vertices of an equilateral triangle. Even 
for such a special case, finding the extrema of the spec- 
tral distribution of photons is unlikely to be carried out 
by hand, and one needs to rely on numerical calculations 
in general. It is shown in the following that the large 
separation limit reduces the optimization problem sig- 
nificantly. Therefore, we shall investigate the relations 
between interference visibility and various known entan- 
glement measures for three-atom systems in this large 
separation limit. The full analysis including resonant 
atom-photon interactions is to be studied in a future pub- 
lication. 

Consider the simple geometry where the atoms are 
pinned at the vertices of an equilateral triangle with 
equal mutual distance r. Take their positions at Xj — 
r(siri$j, — cos$j, 0)/y/3 with the angle dj = 27r(j — l)/3 
(j = 1, 2, 3). As for the two-atom case, the initial atomic 
system is the subspace spanned by \egg), \geg), and \gge) 
while the radiation field starts in its vacuum state |0). 
This subspace of the initial state is referred to as the 
subspace Hp as before. The visibility of the far- field in- 
terference pattern is calculated when a photon is sponta- 
neously released in the mode |fc) with angular frequency 
uj = c\k\ and the atomic system ends in its ground state 
\ggg). Moreover, instead of considering a general initial 
mixed atomic state as was done in the previous sections, 
the following family of pure states, known as the W-like 
states, is investigated 

W) = ci\egg) + c 2 \geg) + c 3 \gge) , (28) 



Cg = 1. Without loss 
j can be set as positive 
C3 > 0. The conditions 



with the normalization c\ + c: 
of generality, the coefficients c 
numbers fulfilling c\ > C2 > 
Cj > arc imposed, otherwise the problem reduces to the 
two- atom case when one of them is equal to zero. Such 
a choice is always possible since any state with arbitrary 
complex coefficients Cj exp(i^j) can always be unitarily 
mapped onto the state (|2"51) by applying the following 
tensor product of local phase gates, 





-i</>i 





e -i02 



1 

e"^ 3 



According to the stochastic local operation and classical 
communication classification of entanglement, the state 
Eq. (|28|) belongs to the W-state class and as such its 
tangle is equal to zero [l8j . However, we comment that, 
even for the pure three-qubit system, there is no generally 
accepted entanglement measure which is able to fully and 
uniquely characterize the amount of entanglement. 



A. Superradiance and subradiance effects for the 
three-atom case 

Similar to the case of two atoms, we briefly state the 
results for the interaction of the three-atom system, ini- 
tially prepared in the subspace Hp, with the scalar radi- 
ation field initially starting in its vacuum state |0). As a 
result of the resonant exchange of photons between the 
atoms, the energy degeneracy of this subspace is split in 
two, one superradiant state and two degenerate subradi- 
ant states. The projected resolvent operator describing 
the evolution of the system in the subspace Hp is given 
as 



G P (z) 



1 



U-) 



z — Ml^ 
1 



1- 



12-) 



1 



z — hfi- 



(2-|, (29) 



where the superradiant and the subradiant eigenstates 
(I = 1,2) are given by 

|+> = (\eggO) + \gegO) + \ggeO}) /V3, (30) 
IM = (q e \eggO)+q 2i \gegO) + \ggeO))/V3, (31) 

where q = exp(2-7ri/3). Taking the origin of energies at 
SHujg and absorbing the Lamb-shift correction into a re- 
definition of the individual atomic angular transition fre- 
quency ujq, the superradiant and subradiant eigenvalues 
fl± = lu± — iT±/2 read 



r 3±l 



r+ = r i ± 



4 

3± 1 



(32) 
(33) 



where / = f(k r) and g = g(k r) are given in (fl4l . As 
a result, the superradiant state |+) decays three times 
faster than the isolated atoms and the degenerate subra- 
diant states \£— ) does not decay in the limit r — > 0. For 
atoms far apart r — > oo, the limits T± — > T and ui± — > loq 
hold, and we recover the case of isolated atoms. 

The initial state \ip)\0) living in the subspace Hp de- 
cays to the atomic ground state \ggg) with a sponta- 
neously emitted photon |fe). In a way similar to the two- 
atom case, the probability of such a process reads 
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P(mo) ->■ \gggk)) 
The weights of the Lorentzians read 

o 3 



Tuj 



( w - w+ )2+(r + /2)2 ( w - w _)2+(r_/2)^ 



£>+ = -c 2 + {c[c + (ci + Cj — 2c)/i + ] cos k-(xi — Xj) + c(a — Cj)ha sinfc- (xt — Xj)} 



(34) 



(35) 



i>j=l 



1 3 

V- = -(1 — 3c 2 ) + {[(cj — c)(cj — c) + c(c, + Cj — 2c)/i_] cosfe-(a;j — Xj) + c(cj — Cj)ho sink-(xi — Xj)} , (36) 



i>j=i 



r 



where c = (ci + c 2 + C3)/3, and the coefficients /i±,o are 
2+/ A . 3±1 



h = 



(3 ff ) 2 + (2 + /)* 
3ff 

(3.g) 2 + (2 + /)2 



1± 



-/ 



(37) 
(38) 



B. Visibility in the large separation limit 

The previous expression for the transition probabil- 
ity (|34[) is quite complex and the large separation limit 
kor ^ 1 is considered to simplify the results and the 
discussion. In this case, the simplified result reads 



P(W)\0) -> \gggk)) 



Tuj 



2(2ir) 2 k$ {uj-ujo) 2 + (r/2) 2 



1 + 2 CiCj cosk-(xi — Xj) 

i>j=l 



(39) 



Indeed, as already discussed in the two-atom case, in the 
large separation limit the radiative decay from state \tp) 
is equivalent to a triple slit experiment where the pho- 
ton could have been released with equal probability from 
any of the states \egg), \geg), and \gge). This is because 
the resonant exchange of photons between the atoms is 
negligible and the degeneracy in the subspace Hp is not 
lifted. The superradiant and subradiant states are de- 
generate in this limit and their radiation properties are 
identical. Note that the coefficient 2ciCj is the concur- 
rence of the reduced state = Trfc^ij |V')(V'I obtained 
when the remaining atom has been traced out. 

To extract the visibility of the fringe pattern, the tran- 
sition probability (|39[) needs to be maximized and min- 
imized over the emission angles of the released photon. 
This is equivalent to finding the maximum and minimum 
values of the function 



1(6 j) = 1 + 2 (c 2 c 3 cos 9i + c 3 ci cos 6*2 + cic 2 cos# 3 ) 



over the angles B. L = k ■ (xj — Xk) 



(40) 



r(cos -dj, siw&j, 0). 



Here (ijk) is a cyclic permutation of 



(123). By definition, these angles are linearly dependent 
as constrained by 9\ + 8 2 + #3 = 0. With this function 
I(9j), the visibility is expressed as 

V(iV))= ^ maJC ~^ min . (41) 

-'max + -'min 

The maximum is easily found when all cosine functions 
reach the maximum by cos 9j — I (j = 1, 2, 3), leading to 

/max = (Cl +C 2 +C 3 ) 2 . (42) 

For example a bright fringe is obtained when the photon 
is emitted perpendicularly to the plane of the equilateral 
triangle made by the three atoms, in which case k-(xi — 

Finding the minimum requires a little bit of caution as 
there are two cases. In the first case, 7 mm = 0, leading 
to full fringe visibility V = 1. As ci > c 2 > C3, this is the 
case if and only if the triangle inequality c\ < c 2 + C3 is 
fulfilled. This is achieved by the angles satisfying 



cos 6j 



2c? 



1 



2ciC 2 C3 

On the other hand, if c\ > c 2 
taincd when cos &\ — 1 and cos 
to 



Cj (3 = 1,2,3). 



(43) 



- C3, the minimum is at- 
->2 = cos 9s = — 1 , loading 



(ci - c 2 - c 3 y 



(44) 



The visibility in this case is less than unity. 

Summarizing the previous results, the visibility reads 



V(|V» = 



1 

2ci(c 2 +c 3 ) 

1 + 2C2C3 



< 1 



(Ci < C 2 + C 3 ) 
(Ci > C 2 + C 3 ) . 



(45) 



We remark that the so-called W-state \W) = (\egg) + 
\g e g) + lff5 e ))/v / 3 achieves V = 1 as naively expected. 
This is easy to understand in the Young slit language as 
the triple slits radiate in this condition on equal footing. 
However it is stressed that there are infinitely many other 
states reaching full fringe visibility. Finally, note that the 
other extreme limiting case c 2 , C3 — > means that only 
one slit is operating and thus one obtains V = in this 
case at least in the large limit. In the next section, these 
results are examined in the light of known entanglement 
measures. 
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C. Relation between interference visibility and 
entanglement measures 

Bearing in mind the results of Sec. Ill, it is natural to 
expect that the previous three-atom interference visibil- 
ity (|45|) is somehow related to the amount of entangle- 
ment present in the initial atomic W-like state ([28)) . Un- 
fortunately, it is seen that W-like states with c\ < C2 + C3 
cannot be discriminated at all since they all lead to V = 1 . 
This raises the question of finding the range of entangle- 
ment measures compatible with a given observed value 
of the interference visibility. Generally speaking, one can 
infer the amount of entanglement in the state with better 
precision when the range of corresponding entanglement 
measure is smaller. Depending upon the protocols used, 
one may need to know the least amount of entanglement 
to make sure that the initial state contains enough en- 
tanglement for the protocols. In these latter situations, 
the lower bounds play a more important role than the 
ranges themselves. 

Our problem is thus to find the maximum and the min- 
imum values of the known entanglement measures com- 
patible with a given value of the interference visibility 
subject to the condition that the initial atomic state is 
the W-like state \ip) given by j28j with c\ > c 2 > c 3 > 0. 
To focus on genuine tripartite entanglement, the con- 
dition C3 > is reminded. When the maximum and 
minimum do not exist, the supremum and infimum are 
calculated, respectively. 

In the following, we study the relations between the in- 
terference visibility and (a) the mixedness of the subsys- 
tem, (b) the geometric measure, (c) the largest bipartite 
negativity, and (d) the three- it. The definitions of these 
entanglement measures are summarized in the Appendix. 
These particular measures are chosen to obtain analyti- 
cal results but other entanglement measures could have 
been studied along the same lines as well. Because these 
calculations are rather tedious, only the final results are 
shown and details will be published elsewhere. 



1. Analytical results 

(a) Mixedness of subsytem: The mixedness of sub- 
system M of the W-like state \tp) is given by 



M 



8 



/ 2 ^ 1 2 a , 2 2 \ 
\ c l c 2 C 2 C 3 c 3 c l) 



(46) 



When ci < C2 +C3, V = 1 and the range of mixedness is 
2/3 < M < 8/9, where the upper bound is attained with 
the W-state and the lower bound is attained by states 
with ci = c 2 + c 3 . 

When ci > c 2 + c 3 (V < 1), the range is 



1 



A 1 



(3 



< M < 



1 + V 2 



(47) 



where v = yl — V 2 is the maximal complementary quan- 
tity to the visibility. The upper bound is obtained for 



C3 — > and this case is excluded. The lower bound is ob- 
tained for c 2 = c 3 . Interestingly, these bounds converge 
to the same value 2/3 when one takes the limit V — >• 1 
from below. 

(b) Geometric measure: The geometric measure E g 
for three qubits has been studied for the above W-like 
states (|28|) and has been obtained analytically as a func- 
tion of the coefficients Cj If c\ < c 2 + C3 (which 
implies c\ < y/2/2), then c\ < c 2 + c 3 and one can con- 
struct an acute triangle with the coefficients Cj being the 
lengths of its edges. In this case 



E. 



g ~ 
R = 



1 -4ir 



ClC 2 C 3 



4\/co(co - ci)(c - c 2 )(c - c 3 ) 



(48) 
(49) 



where R is the circumradius of the triangle formed by 
the Cj and cq = (c% + c 2 + c^)/2. In the other case [i.e., 
c i > c 2 + c i ( c i > \/2/2)] one finds 



Eg = 1 - 4 . 



(50) 



When c\ < c 2 + C3, the range of the geometric measure 
is obtained as 1/3 < E g < 5/9, where the upper bound 
is found for the W- stat e and the lower bound is found for 
states with c\ = \/2/3. 

On the other hand, when ci > c 2 + C3, the geometric 
measure is bounded by 



1 



v _ 1 — v 

— < Eg < - 



(51) 



where the upper bound is obtained for states with C3 — > 
and the lower bound is obtained for states where c 2 = C3. 
These upper and lower bounds asymptotically reach the 
values 1/2 and 1/3, respectively, asV->l from below. 

(c) Largest bipartite negativity: The bipartite neg- 
ativity of the W-like states with respect to the jth qubit 

is given by Afj = 2cj*\ — c 2 . The maximum value 



achieved over all possible bipartite partitions is of interest 
and thus we consider the quantity 



AC 



:{Ai,A/2,A3} 



(52) 



When ci < c 2 + c 3 , the range of the maximum bipartite 
negativity is 2 v / 2/3 < Af max < 1- The lower bound is 
found for states with c\ = y/2/3 or c\ = \/l/3 while the 
upper bound is found for the W-state. 

When ci > c 2 + C3, the largest bipartite negativity is 
bounded by 



1 + 3v 



< ■N'ma.x < V . 



(53) 



where the upper bound is found for states with C3 — > 
and the lower bound is found for states with c 2 = c 3 . The 
upper bound reaches 1 and the lower bound converges to 
v8/3 (= 0.9428 . . . ) as V approaches 1 from below. 
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(c) (d) 

FIG. 2: The ranges of entanglement measures studied in this paper as a function of the visibility of the far-field interference 
pattern. The atoms are pinned at the vertices of an equilateral triangle, and their mutual distance is large, (a) Mixedness of 
subsystem, (b) geometric measure, (c) bipartite negativity, and (d) three-7r, respectively. In each graph, the range for V < 1 is 
shown by two curves filled with a gray area, and the thick vertical line indicates the range at V = 1. The asymptotic values for 
V —¥l from below are shown by the dashed-dotted lines. 



(d) Three- ir: All entanglement measures studied pre- 
viously do not vanish when C3 — > 0. In other words, they 
are insensitive to genuine tripartite entanglement, but 
are rather related to entanglement stemming both from 
bipartite and tripartite entanglements. This is seen from 
their definitions that all three entanglement measures do 
not vanish when the problem reduces to the two-atom 
case (i.e., c\ > c-i > C3 = 0). 

To examine the relation to genuine tripartite entangle- 
ment measure, the so-called three-7r is examined, which 
has a similar origin as the tangle 0, l2lj |. This is ex- 
pressed in terms of the Cj as 

^ = ^cj(/iT^^4-l) . (54) 

As is easily checked, Af v = holds when C3 — > show- 
ing that this entanglement measure is only sensitive to 
genuine tripartite entanglement. 

When ci < C2 + C3 the range of the three- ir is < 
AC < 4(a/5- l)/9, where the upper bound is attained 
by the W-state and the infimum is obtained for states 
with C3 — > 0. When c\ > C2 + C3, the three- n is bounded 
by 

< AC < AC iax , (55) 

where the infimum is again obtained for states with C3 — ¥ 
0. The upper bound is found for the state with ci — C3 = 



Vvl — w )/2(3 + v) and reads 

= iTTT^lW 1 + ^5 + 6^ + 5^ 

3(3 + vy 

+ (1 - v)y/(l - v)(17 + 5v) - (9 + Uv + 9v 2 )] . (56) 

The upper bound converges to the value 2(4-\/5 + VTf — 
9)/27 (= 0.3012 . . .) as V -> 1 from below. 

The results are plotted in Fig. 2, (a) mixedness of sub- 
system, (b) geometric measure, (c) bipartite negativity, 
and (d) three-7r, respectively. In each graph, the range 
for V < 1 is shown by two curves filled with a gray area, 
and the thick vertical line indicates the range at V = 1. 
The asymptotic values for V — > 1 from below are shown 
by the dashed-dotted lines. 



2. Summary of the results 

Let us summarize the obtained relations between the 
fringe visibility and the four entanglement measures and 
discuss some consequences of them. 

First, the mixedness of the subsystem and the geo- 
metric measure exhibit similar behaviors showing jumps 
in the ranges when V — > 1 from below. These jumps 
also exist in the three- 7r case, and result in detecting the 
amount of entanglement around the full visibility with 
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less accuracy. The appearance of these jumps is partly 
due to the fact that interference visibility defined by ^ 
is a quantity closely related to a bipartite system whereas 
the mixedness of the subsystem and geometric measure 
are related to both bipartite and tripartite entanglement, 
and the three- 7r is only related to tripartite entanglement. 



Second, Fig. 2(c) shows that interference visibility can 



reasonably be used to detect the amount of bipartite en- 
tanglement in the W-like state for all values of interfer- 
ence visibility. Remarkably, the upper bound is exactly 
equal to the largest amount of bipartite negativity with 
respect to all possible bipartite separations and does not 
show any jumps when the visibility approaches its max- 
imal value V = 1 from below. 



Last, Fig. 2(d) clearly indicates that the interference 
visibility is a poor witness of the genuine tripartite entan- 
glement initially present in the atomic system. In par- 
ticular, getting a far-field interference pattern with full 
visibility V = 1 does not tell anything at all about the 
least amount of tripartite entanglement initially present. 
This result is not really surprising since ([59")) shows that 
the interference pattern directly depends on the prod- 
ucts CiCj. As mentioned, these products identify with 
the concurrence of the two-atom reduced states obtained 
by tracing out one atom. It is thus obvious that the 
interference pattern perfectly encodes some bipartite en- 
tanglement properties of the atomic system, but is less 
sensitive to genuine tripartite entanglement properties. 



D. State tomography for the W-like states 

In this last section, a state tomography is briefly dis- 
cussed for the W-like states ([28]) from information con- 
tained in the far-field interference pattern by observing 
spontaneously emitted photons. The coefficients Cj can 
be extracted simply by taking a few points of the inter- 
ference pattern l|39p. Because of the normalization condi- 
tion of the state, one needs at least three different points 
to reconstruct the state. 

When the W-like initial state has complex coefficients 
Cj — > Cj exp i<j)j , the full state tomography also needs to 
find the phases cf>j. Since the global phase of the state 
is irrelevant c\ can set to be real. The remaining phases 
4>2 and </>3 are incorporated in the interference pattern 
through the replacements 9j — > 8j\ Q\ = Q\ + <j>2 — 4>3, 
&2 = 82 — 03 ! #3 = 83 + <\>2 ■ This shows that these phases 
(j>j are merely responsible for a shift of the interference 
pattern. Knowing the position of the fringes, one can ex- 
tract 4>2 and 4>3, thus completing the state tomography. 
We comment that a more complicated analysis is needed 
when the initial atomic state is mixed. The full tomogra- 
phy may well be possible when the atoms are separated 
by a distance on the order of optical atomic transition 
wavelength. 



V. CONCLUSION AND OUTLOOK 

In this paper, a system of two-level atoms has been con- 
sidered when initially prepared in the subspace spanned 
by the first excited atomic states without photons and 
releasing a spontaneous photon in the course of time. 
Some detailed analysis has been given about what kind 
of information the observation of the far- field interference 
pattern can help reveal on the amount of entanglement 
present in the initial atomic state. 

In the two-atom case, a previous known result is con- 
firmed showing that the interference visibility indeed con- 
verges to the concurrence of the initial atomic state pro- 
vided the two atomic qubits are far apart. For finite 
relative distances between the qubits, the deviation of 
the visibility from the concurrence is calculated numer- 
ically. It is found that this deviation is small whenever 
the qubits are separated by a distance r which is a mul- 
tiple integer of half the wavelength Ao associated to the 
optical atomic transition. A complete state tomography 
about the state is shown to be possible provided the two 
atoms are sufficiently close to each other. 

In the three-atom case, the family of W-like states has 
been considered and the relations between the interfer- 
ence visibility and several known entanglement measures 
have been analyzed. Bounds have been derived analyti- 
cally for the amount of entanglement present in the initial 
state when the interference visibility is given. Our results 
show that the interference visibility is closely related to 
the bipartite entanglement properties of the atoms. In 
particular, it is shown that the interference visibility is 
directly related to the maximal bipartite negativity. 

Two possible extensions of this work are to be men- 
tioned. First, the possibility of detecting genuine multi- 
partite entanglement by defining a suitable interference 
visibility for multipath interferometers. Second, the anal- 
ysis of the wave-particle duality aspects in the three-atom 
case when the atoms have an internal structure allowing 
for the storage of which-path information. These issues 
shall be investigated in the future. 
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Appendix A: Entanglement measures 

Definitions of the entanglement measures used in this 
paper are listed. All these measures take on values be- 
tween and 1 as a convention and are defined for two- 
qubit and three-qubit systems living in the Hilbert space 
(C 2 )®" (n = 2,3) of the forms @7|) and (J3SJ). For more 
precise definitions of these entanglement measures and 
for an account of their properties, we refer to Refs. [22|,[23[ 
and references therein. 



1. Concurrence 

The concurrence of a general two-qubit state p is de- 
fined by 

C{p) = max{0, VA7 - V^i - \/h - . (Al) 

where Ai > A 2 > A3 > A4 are the eigenvalues of yfp {a^® 

cr y 2 ^)p* {<j\P ® &]p)-i/P- The y component of the Pauli 

spin operator is written as <Ty — — i\e)(g\ + i\g)(e\, and 
p* denotes the complex conjugate of the state p. 

2. Negativity 

The negativity of a bipartite system, described by a 
state p, is defined as twice the sum of the negative eigen- 
values of the state p Tj obtained by partial transpose with 
respect to one subsystem. This can be formally expressed 
as 

A^) = ||P T1 ||-1 = IIP T2 H-1, (A2) 

where the symbol Tj denotes the partial transpose with 
respect to the jth subsystem and ||X|| is the trace class 
norm of the operator X. By convention the negativity of 
the two-qubit Bell states is set to 1. 

For multipartite system, the negativity is defined by 
the bipartite partitions of the system. For example, 
■Afj(p) = \\p Tj \\ — 1 can be defined to quantify the en- 
tanglement between the jth qubit and the rest of the 
system. 

3. Mixedness of subsystem 

The mixedness of the subsystem as an entanglement 
measure in a pure state p = is defined by the 



arithmetical mean of the mixedness associated to each 
one-qubit reduced states obtained from p. For the three- 
atom case for example, it reads 

M(p) = 1 -{S{ P1 ) + S( P2 ) + S( P3 )) , (A3) 

S(p j ) = 2(l-Trp 2 j ), (A4) 

where p\ = Tr2,3 p, and so on. In several papers, this en- 
tanglement measure coincides with the definition of visi- 
bility (e.g., Ref. 0,0,1 El). 



4. Geometric measure 

The geometric measure for an n-qubit system de- 
scribed by a pure state \ip) is defined by the distance 
between \ip) and the closest product state 

E g (p)= min ||^)-|0)| 2 , (A5) 

|0)£Pro 

where Pro represents the space of product states within 
the Hilbert space C® n . 



5. Three-7r 

Several authors have introduced and studied this mea- 
sure by following the idea which led to the tangle in 
Ref. [20|. The most convenient construction appeared 
in Ref. [2l| with the definition 

K(P) = ^(TTl +7T2+7T3), (A6) 

TTj = A/» + N 2 { Pj ) j^^iPj) , (A7) 
j'=i 

and Af(pj ) is the negativity of the reduced state over the 
jth qubit [i.e., pj — Trj(p). This can be written more 
compactly as 
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